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ABSTRACT 

It is shown that squeezed states of primordial gravitational waves are inevitably produced in the 
course of cosmological evolution. The theory of squeezed gravitons is very similar to the theory of 
squeezed light. Squeeze parameters and statistical properties of the expected relic gravity-wave radia- 
tion are described. 

INTRODUCTION 

Squeezed quantum states of light have been successfully generated and detected under laboratory 
conditions. It is known how much skill and effort by our experimentalist colleagues it requires to 
achieve even a modest amount of squeezing, that is, to obtain the squeeze parameter r of order of 1. 
The main purpose of my talk is to show that, in the cosmos, squeezed quantum states of gravitational 
waves are produced inevitably and with a much much larger amount of squeezing, simply as a result 
of expansion of the Universe. 

In the context of gravity-wave research, the notion of squeezed quantum states has been often 
referred to. However, what was always meant was the squeezing of a quantized vibrational mode of a 
detecting device that could be implemented for a better detection of a classical gravitational wave. For 
instance, it was shown (Ref. 1) that the performance of a laser interferometer gravity-wave detector 
can be improved by using squeezed light. In another paper (Ref. 2) it was argued that any detector- 
oscillator can be specially “prepared” in a squeezed state and used for gravity-wave detection during 
some interval of time before the thermal noise destroys squeezing and degrades the detector’s sensitivity. 

However, it is the squeezing of the gravitational waves themselves that will be discussed in my 
talk today. I will show that the production of squeezed relic gravitational waves is an inescapable 
consequence of the variability of cosmological gravitational field and the existence of zero-point quantum 
fluctuations. 

The mathematical theories of relic graviton production and squeezing of light are very similar. To 
make this similarity especially transparent, I will begin by presenting Einstein’s general relativity in the 
form of a traditional field theory, such as the theory of classical electromagnetic fields. Those of you w r ho 
may feel uncomfortable, or even intimidated, with the notion of curved space-time, will, perhaps, find 
it easier to deal with the concept of a gravitational field given in the usual flat Minkowski space-time. 
(More details about this “fi eld- theoretical” formulation of general relativity are presented in Ref. 3; it 
is important to emphasize that we are dealing with a different mathematical formulation of general- 
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relativity, not with some alternate physical theory, see Ref. 4.) This approach leads to manifestly 
nonlinear field equations. In contrast to quantum optics based on the laws of linear electrodynamics, 
in the case of a gravitational field one does not need any material medium in order to couple a “pump” 
field with the “signal” waves; for gravity this is achieved automatically due to the nonlinearity of the 
gravitational field itself. As is often done, we will present the total gravitational field in the form 
of an approximate sum of a large “classical” contribution and a small quantized perturbation. This 
approach will be applied to the cosmological gravitational field of the expanding Universe acting upon 
zero-point quantum fluctuations of the gravitational waves. I hate to call the variable gravitational field 
of the expanding Universe -the most grandiose and magnificent phenomenon we are aware of - just a 
“pump” field, but, technically speaking, it plays precisely this role. As a result, the initial vacuum state 
of gravitational waves will evolve into a strongly squeezed vacuum state with very specific statistical 
properties. 

We will discuss the expected characteristics of the relic gravity-wave background radiation and the 
problem of its detection. 

FIELD-THEORETICAL APPROACH TO GENERAL RELATIVITY 


A gravitational field is fully described by a symmetric second-rank tensor (note that the gravita- 
tional field variables have just one extra index as compared with the electromagnetic 4-vector potential 
A M , not a big difference!). For writing down the Lagrangian of the gravitational field it is convenient 
to use also an additional set of variables: the tensor field P a ^, symmetric with respect to the last two 
indices. However, is not a new physical field but rather a combination of the first derivatives of 
as follows from the field equations. 


Gravitational field potentials h fil/ (x y y J z^t) are mathematically treated as components of one of 
many physical fields immersed in the ordinary Minkowski space-time: 

da 2 — c 2 dt 2 - dx 2 - dy 2 — dz 2 . (1) 

The metric tensor of Minkowski space-time will be denoted by 7 Mt/ . With respect to this metric tensor 
all covariant differentiations (denoted by the symbol and lowerings or raisings of indices are to be 
performed. In the Lorentzian coordinates, like the ones implied in eq. (1) and which we will be using 
in practical calculations below, 7 MI/ acquires especially simple values: 700 = 1, 7n = 722 = 733 = — 1, 
with the rest of 7^ being equal to zero. 


The gravitational part of the total action S = S 9 + S m is 


S 3 = — f d 4 x L 9 , k — 

2 CK J 

where the gravitational Lagrangian L 9 has the form 


8 ttG 


L 3 = (-7 ) 1/2 [h^, Q p a „ - {hr + ri (/>“ 
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The nongravitational sources and fields and their interaction with the gravitational field are described 

by 

S m = ~ c j d 4 x L m . 
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The energy-momentum tensor t of the gravitational field itself and the energy-momentum tensor 
of the nongravitational matter and fields interacting with gravity are defined in the usual manner: 


1 6LP 2 6L m 

Ki "" ~ "7=7 ~fr ' Tilv ~ 7=7 tv* ■ 

The precise expression for is as follows: 

— P a ^pP^ — ~7, ^1° P P° * trpP^pa ” g ^ 

where P» = = Q a M „ ;a and 

<?%„ = \p ° a0 [i^ni + 7 ™ ( 7 „x + («: t™ + - « (<:* f + • 


By varying the action S with respect to and one can derive the following gravitational 

field equations: 

\ h-ir/J 

(3) 


167tG 


hftu ’ ,a "P l P u h ,a,0 h v,a,n ^ ~ ^4 (^M*' P 7 V*') 


where the comma means an ordinary derivative (, a = d/dx a ) and 7„„ is assumed to be in the simplest 
form corresponding to eq. (1). For reference purposes we will also write down the relationship between 
the first derivatives of h pu and P a M „: 


-hr p -h^r 0 Ia0, fi +(r a + v a ) p v qp Hi vo + hl,Q ) p ^ P -\Pa [(r a + h^s; + ( 7 *'°+ hr)%\ = 0 . 


The theory possesses a gauge freedom quite similar to the gauge freedom of classical electrodynam- 
ics. One can apply the gauge transformations to the gravitational variables h pi/ and matter variables 
without changing the field equations. At the expense of the gauge freedom one can impose some gauge 
conditions which are normally used for diminishing the number of variables and simplifying the field 
equations. 

The transition to the usual “geometrical” formulation of general relativity is established by intro- 
ducing the new functions g pi/ according to the rule 

7=5 <T = 7=7 (V' + *'"') (4) 

and by identifying the g pv with the metric tensor of the curved space-time: ds 2 = g Q0 dx a dx 0 . 

GRAVITATIONAL FIELD OF THE EXPANDING UNIVERSE 

Let us apply the developed formalism for description of the gravitational field of the homogeneous 
isotropic Universe. From our new point of view this is just a specific gravitational field h^(t,x,y,z) 
given in Minkowski space-time (1). Let us take the nonvanishing gravitational potentials in the form 

hoo = a 3 (<) - 1 , hn = ^22 = /»33 = 1 — a (t) (5) 
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where a(t) is, as yet unspecified, function of time. One can calculate the gravitational energy- 
momentum tensor t ul/ , eq. (2), and find that the nonvanishing components of are 


Ktoo = -^(a 2 -l)-3a 2 
2a 

k<h = k* 22 = = -^-(3a 3 - a 2 + a - 3) - a 2 (3a - 1) (6) 


where the dot means the time derivative (for simplicity we choose the units where the velocity of light 
c = 1). 


The nongravitational sources are assumed to be in the form of hydrodynamical matter with the 
Lagrangian 

L m = \y/~9 [e + 3p - (e + p)g lll/ tt" u v ] 

where g is defined by eq. (4). One can find the nonvanishing components of the energy-momentum 
tensor 


r oo 


Til = r 22 = T33 


c-t-^(a 2 - 1 ){(~p) 

P ~\( a2 ~ l )( ( - P ) ■ 


(7) 


By substituting expressions (5), (6), and (7) into the field equations (3) one can derive equations 
governing the function a(t) and, hence, the gravitational field (5): 


a 

a 


4ttG 


(c + 3p), 



8 t vG 
€ 


(In “geometrical” language, these are, of course, the Einstein equations for a spatially flat cosmological 
model: ds 2 = dt 2 — a 2 (t) ( dx 2 + dy 2 -f dz 2 ).) By specifying the relationship between e and p (“the 
equation of state”) one can solve these equations and find a concrete function a(t). 


AMPLIFICATION OF GRAVITATIONAL WAVES 


Gravitational field (5) is just a main term of a more complicated and realistic cosmological gravi- 
tational field which includes the gravity-wave perturbations. Let us present the total field in the 
form 

= 1$ + *£?(*» x > Vi z ) - ( 8 ) 

where hf$ is given by eq. (5). At the expense of gauge freedom one can always satisfy the conditions 
_ q Moreover, in the case of gravitational waves, the perturbations of e, p , and u a are equal to 
zero and one can, in addition, satisfy the requirements = 0, h$ 7 MI/ = 0, so that one is left with 
only two independent polarization components (designated by a = 1, 2 ) of h$. For a wave with the 
wave vector n one can write down the nonzero components of the field: 

h\] ] (t, x ,y,z) = (/£(*)«’“ + /C(0 e_,nx )po (9) 

where the constant polarization tensors p? fulfill the conditions p“ fc n fc = 0, p“i = 0. 
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Now one should substitute (8) into the field equations (3) and linearize them with respect to 
It is clear that the left-hand side of eq. (3) is simply the usual D'Alembert differential operator applied 
to At the same time, the right-hand side of eq. (3) will contain the products of and h$ since 
all the nonlinearities are collected there. 

For a given perturbation with the wave vector n and for each of the two polarization components, 
the field equations reduce to a single equation for the time-dependent function ^i(t) (indices n and a 
are omitted): 

fi + n 2 fi = 

where n 2 = n* 2 + n 2,2 + n 32 . If there is no “pump” field (5), that is a(t) = 1, the right-hand side of eq. 
(10) vanishes. It is worth noting that in the gravitational case, in contrast to electrodynamics, there 
is not any dimensional coupling constant between the “pump” and “wave” fields, the strength of the 
coupling is regulated by the rate of variability of the “pump” field. 

It is convenient to introduce a new time coordinate 7] related to t by drj = a(t)~ l dt , and to denote 
the r \- time derivative by prime. Equation (10) gets an especially simple form (Ref. 5): 



which makes it possible to treat the problem as a problem for a parametrically excited oscillator. The 
notion of squeezing appears quite naturally. 

SQUEEZED VACUUM STATES OF RELIC GRAVITONS 

In some cosmologically interesting and realistic situations, the function a" fa goes asymptotically 
to zero for r] — ► -oo and rj — ► + oo. In the asymptotic regions rj — ► — oo and r) — ► + oo, solutions to eq. 
(11) are very simple: fi(rj) ~ e itnf7 The general complex solution to eq. (11) can be presented in the 
form 

Kv) = «£(»?) + * + C0?) ( 12 ) 

where £(r)) and C(v) are complex-conjugated normalized base functions. The same general solution 
can be decomposed over other base functions x(t)) and 

Kv) = c x(v) + <Px*(v)- ( 13 ) 

One can choose the base functions in such a way that 

t(v) -* ~ir= e ~' nT ' for -oo 

V2 n 

and 

xin) — -m =e~' nv for tj +oo . 

V2n 

Since (12) and (13) describe the same solution, their coefficients are related: 

a = uc + vd+ , 6 + = v*c + u*d^ , (14) 

where |u| 2 — \v\ 2 — 1. 


n 2 (a 2 — 1 ) d 


+ “+ ( - 
a \a 
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For a quantized field, the coefficients a, 6 + , c, have the meaning of the creation and annihilation 
operators and the relations (14) are called the Bogoliubov transformations. 

Complex numbers u,v can be parameterized by the three real numbers r, 0, </?, r > 0: 

u — e~ %d chr , v = — e x ^+ 2 ^shr . (15) 

The transformations (14) can also be presented in the form: 

a = R+S+cSR, 6+ = R+S+d+SR , (16) 

where £(r, <p) and R(0) are unitary operators: 

S(r,<p) = exp [r(e~ 2t ^cd — e 2tv c+d + )] 

R(0) = exp[— i0c + c - i9d+d\ . 

In the theory of squeezed quantum states, the operator S(r y <p) is called the two-mode squeeze operator, 
the operator R(9) is the two-mode rotation operator, r is the squeeze parameter and is the squeeze 
angle (see, for instance, Ref. 6). 

As a result of evolution, the two- mode vacuum state |0,0) transforms into a two-mode squeezed 
vacuum state: 

\SS ) 2 = 5(r, V )|0,0) . 


The two modes under discussion are two waves with the same frequency but propagating in opposite 
directions. In the field of quantum optics, one is usually interested in the temporal fluctuations of the 
light field, so that the spatial distribution of the field is not always important. However, in cosmology, 
we need to know the complete space-time distribution of the gravity-wave field. For this aim we augment 
the time-dependent functions fi(rj) with the spatial functions t/(x). For every n-mode contribution and 
a given polarization component, one has: 

h n = Mn, l Un,\ + /V2^n,2 * 

One can work, for example, with complex functions Un , l = K e tnx , U n ,2 — K e“* nx (and complex 
conjugated /z: p n , 2 = Mn,i) or rea l functions t/ n , 1 = K cosnx, f / n ,2 = K sin nx (and real ji), 
where K is a normalization constant. Classically, this corresponds to the decomposition of the field 
over traveling or standing waves. For the field operators, one writes 

h n = « + t> + C) A'e ,nx + (a+r + K e~' nx , 


in the first case, and 

h n = (tnt + b+c) V 2 K cos nx + fa + b+C) V2K sin nx , 
in the second case. Transition between the two descriptions is fulfilled by the transformation 


CL -f- b 


= 6 , 


i(a - b) 

V2 


= 62 • 


(17) 


(18) 


In terms of the theory of squeezed states, it means that, in the second case, one will be dealing with a 
pair of one-mode squeezed states instead of a single two-mode state. Indeed, under the transformation 
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(18) the two-mode squeeze operator factorizes into a product of the two one- mode squeeze operators 
, so that instead of (16) one will have 

bi = RfStaSiRi 

where 5i(r, <p) is the one-mode squeeze operator 

Si(r,<p) = exp ^( e- 2,ifi c \ - e 2 "*c\ 2 

and R\{9) is the one- mode rotation operator 

Ri(0) = exp(— iffc+ci) . 


The operator b 2 transforms in precisely the same manner and is associated with the other squeezed 
mode. 

As a result of evolution, the one-mode vacuum state transforms into a one-mode squeezed state. 

In what follows, we will base our analysis on the representation (17) and on the one-mode squeezed 
states. The annihilation and creation operators will be denoted by b and 6 + . 

The classical equations of motion (11) can be derived from the Hamiltonian 


H = \ p 2 + — ( up + p/i) + n 2 p : 

4 A 


where p is the momentum, canonically conjugated to the coordinate /i, p — p f — In a standard 
manner, by introducing the annihilation and creation operators 6, 




one can present the Hamiltonian in the form 

H = nb + b + a(Ti)b + ' + a*(ji)b 2 (19) 

where the coupling function <t(t 7) is <7(77) = ia'/2a, and the Planck constant h = 1. Note that the 
Hamiltonian (19) belongs to the class of Hamiltonians that characterize a number of physical processes 
(Ref. 7). However, in most of them the function <r(t) has a specific form a (t) = ae~ 2,UJt where a is a 
constant, albeit in our case a(rj) is a more general function of time. 

The Heisenberg equations of motion, following from this Hamiltonian, have the form 

. db .a' . ,db + . .a' 

1— = nb + 1 — £r — j— — = nfr — 1 — 6 . 

v-i 1 Win rt 


Their solution is 


b(rj) = u(r})bo + v(r))b£ , b + (r}) = u*(i ])b£ + v m (r))bo , 
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where &o, &o are i ni ti aJ values of 6(r/), 6 + (r/) (Schrodinger operators) and the complex functions 

u(t]) and v(r]) satisfy the equations 

iu f = nu + i—v* , —iv** = nv* - i—u , u(0) = 1, u(0) = 0 . (20) 

a a 

It follows from these equations that u + v* satisfies the equation identical to eq. (11): 

( u + v*)" + (n 2 - — )(u + v*) = 0 . 
a 

As for the function u — v*, it can be found from the relation 


-in(u - t;*) = (u + v*)* (u + v*) . 


By substituting eq. (15) into eq. (20) one can find equations for the time-dependent parameters r(7/), 


r 


/ 


o' = 




t 


cos 2(f 

a 

a! , 

n sin2<y? thr 

a 

o! 

—n sin 2 <p cth2r . 

a 


( 21 ) 


Solutions to these equations determine the precise form of evolution of the initial vacuum state into 
a one-mode squeezed vacuum state. Statistical properties of the final state depend on the numerical 
values of r and <p in a well-known way (see, for instance, Refs. 6,7). 


A possible way of calculating r(r?) and <p(rj) (Ref. 8) for a given gravitational field a(r}) relies on 
the observation that the complex function B{r ) ), where 


fu — v*\* _ chr + e 2 "*shr 2 

Vti + v*/ ~ chr - e 2v *shr n ^ ’ 


( 22 ) 


satisfies the equation 

73 2 n f 

B' = i 2-B -2iB 2 

2 a 

with the solution 

OtM 

2 /i/a 

where /i(r/) obeys eq. (11). The properly chosen solutions to eq. (11) define B(t /) and allow one to find 
r(r)) and ip(r /) from eq. (22). The meaning of the function B(r]) is that it determines the Gaussian wave 
function ^(/i,7/) ~ exp(~ B(r])}i 2 ) which is a solution to the Schrodinger equation in the coordinate 
representation. 



The parameters r(7/), (p(r /) can be calculated (Ref. 8) for a cosmological model which includes 
three sequential stages of expansion: inflationary (a(£) ~ e Hot ), radiation-dominated (a(J) ~ t}I 2 ) and 
matter-dominated (a(J) ~ t 2 / 3 ). It can be shown that the present-day values of the squeeze parameter 
r range from r % 1, for frequencies u « 10 8 Hz, up to r % 120, for frequencies v & 10~ 18 Hz. In the 
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frequency interval v « 10 _1 - 1CT 3 Hz, accessible for the planned Laser Interferometer Gravitational- 
Wave Observatory in Space (Ref. 9), the squeeze parameter r reaches large values of order 40-50. 

As for the parameter <f >(t/ ), it can be shown to have the form <p « —nrj + <p o, where <po is a constant. 
This behaviour can already be envisaged from eq. (21) for y/, since, in the asymptotic region rj — ► +oo, 
one has |^-| < n and cth2r « 1. 

RELIC GRAVITONS: A STOCHASTIC COLLECTION OF STANDING WAVES 

As we see, the cosnx and sin nx modes in the representation (17), evolve into a strongly squeezed 
vacuum state. The mean number of quanta (N) and its variance ((AJV) 2 ) are determined by the 
squeeze parameter r: 

(N) = sh 2 r, ((AN) 2 ) = ^sh 2 2r . 

The mean values of ft and p are equal to zero, but their variances do not vanish: 

((Ap) 2 ) = —(ch2r - sh2r cos 2<p), ((Ap) 2 ) = ^( ch2r + sh2r cos 2<p). 

2 Ti 2 


In order to relate the rigorous quantum-mechanical treatment, described above, with the notions 
of random classical waves, one can use the Wigner function formalism. It allows one to derive the 
distributions of the random variables A and <f> entering the classical expression for p: 

p = Asin(— nr) + <f>) . (23) 


It can be shown (see Ref. 8) that, in the limit of large r, the Gaussian distribution for <t> is very 
narrow, like a ^-function. It is concentrated near the values 

<f> = <po + wl (24) 

where <po is a constant, the same for all unit vectors n/n, and l = 0, ±1, — 

A similar conclusion can be reached in a simpler, though perhaps less rigorous way. Let us consider 
the ratio ((Ap) 2 ) / ((A p) 2 ). For large r, this number is approximated by 

((Ap) 2 )/ ((Ap) 2 ) « tg 2 (p . 

For a classical expression (23), this ratio corresponds to the number n~ 2 tg 2 (-nr) + <f>). From their 
comparison, and taking into account the fact that ip « -nr/+ cpo, one can derive eq. (24). The very 
small variance of the phase, A<f>, is, of course, consistent with the large variance of the number of 
quanta, ((AiV) 2 ). 

The negligibly small variance of the phase distribution leads to an important result: every pair 
of cosnx, sin nx modes forms together a standing wave. Indeed, let us consider a given n. The 
corresponding terms, contributing to the total wave-field h(r),x,y,z), can be written in the general 
form: 

h n = A\ sm( — nr} + <?!>i)cos nx + A 2 sin (-nr} + <fo)sin nx . (25) 
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The amplitudes A\ and A 2 are taken from a broad Gaussian distribution and are, in general, different. 
However, the phases <f>\ and <f > 2 are taken from a very narrow Gaussian distribution and are essentially 
fixed and equal up to ±tt. Because of that, expression (25) can be written as a product of a function 
of time and a function of spatial coordinates: 

h n = ± sin {-nr) + V ? o)(^i cos nx ± M sin nx) . (26) 

In other words, expression (26) describes a standing wave. A characteristic feature of a standing wave 
pattern is that the field vanishes all over the space at every half of the period. The randomness of the 
wave-field is displayed in its spatial functions A\ cos nx ± A 2 sin nx. This is why we say that the relic 
gravitational waves are present now in the cosmos in the form of a stochastic collection of standing 
waves. 

The total field h(rj,x) is obtained by summing over all n-mode contributions (26). Of course, the 
total field loses the property of vanishing at some moments of time, because the various sin {—nr} + 
y?o) factors have different arguments. However, the difference in the arguments is not random, but 
deterministic. For instance, if at some moment of time 77 = rjo the component /i n (r/, x) vanishes, the 
same will be true for all other components h m (r), x), where m = n(l + kji) and kfl is an arbitrary 
rational number. Hopefully, this property can somehow be used in a specific strategy of observational 
discrimination of relic gravitational waves from stochastic gravitational waves of a different origin. I 
think that the inevitable “squeezing” of relic gravitational waves (and other primordial fluctuations 
of quantum-mechanical origin) can manifest itself in a variety of circumstances, not all of which are 
foreseeable at the moment. 
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